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The interplay of topological constraints and Coulomb interactions in static and dynamic prop- 
erties of charged polymers is investigated by numerical simulations and scaling arguments. In the 
absence of screening, the long-range interaction localizes irreducible topological constraints into tight 
molecular knots, while composite constraints are factored and separated. Even when the forces are 
screened, tight knots may survive as local (or even global) equilibria, as long as the overall rigidity 
of the polymer is dominated by the Coulomb interactions. As entanglements involving tight knots 
are not easy to eliminate, their presence greatly influences the relaxation times of the system. In 
particular, we find that tight knots in open polymers are removed by diffusion along the chain, 
rather than by opening up. The knot diffusion coeflficient actually decreases with its charge density, 
and for highly charged polymers the knot's position appears frozen. 

PACS numbers: 02.10.Kn 82.35.Rs 87.15.-v 36.20.Ey 05.40.-a 



I. INTRODUCTION 

A polymer chain can be easily deformed, but since 
it cannot cross itself, it is subject to topological con- 
straints. These constraints can be temporary, such as 
entanglements between linear polymers, or permanent if 
the chains are closed (ring polymers) or cross-linked. Un- 
derstanding the influence of topological entanglements 
on static and dynamic properties of polymers is a long- 
standing issue [|], ^, which has recently found renewed 
interest in the context of knotted biopolymers. DNA in 
the cell can change its topology by the topoisomerase 
enzymes that pass one strand through another, in the 
process either creating or removing knots Synthetic 
RNA trefoil knots have been used to prove the existence 
of a similar (previously unknown) topology changing en- 
zyme Q. There is also much interest in developing ar- 
tificial biopolymers, for example as molecular building 
blocks or for DNA-based computing, and in this quest 
complex knots and links have been created in both single 
and double stranded DNA Tight knots have been tied 
in single molecule experiments on both DNA and actin 
filaments using optical tweezers 

Several theoretical approaches have addressed the in- 
fiuence of topological constraints in polymer networks 
and solutions. In particular, the tube model Q in which 
the constraints are replaced by a hard confining tube, 
is quite successful in predicting relaxation dynamics of 
polymeric solutions. In a complementary approach, topo- 
logical constraints are described in terms of localized en- 
tanglements or knots, that perform collective motions 
along the polymers 0. Single molecule experiments are 
now able to probe polymers of specified topology, and to 



examine the infiuence of knot complexity on basic physi- 
cal properties such as the radius of gyration Rg. A simple 
scaling picture ||] suggests that Rg is reduced as a power 
of the knot complexity, measured by the minimal num- 
ber of crossings in a projection. Indeed, a Flory mean 
field theory of knotted ring polymers ^ incorporating 
this knot invariant predicts various scaling dependences 
on knot complexity. A topological localization effect is 
also suggested, in which knots segregate in a single rela- 
tively compact domain while the rest of the polymer ring 
expels all the entanglements and swells freely. Recent 
Monte Carlo simulations in Refs. Q support 

the idea that entropic factors localize topological con- 
straints. This is bolstered by analytical arguments on 
slip-linked polymers and experiments on vibrated 
granular chains [ |l6[ . 

Many biopolymers are highly charged. The effect of 
electrostatics on knotting probability of double stranded 
DNA has been studied in the case where the screening 
length is smaller than the persistence length of the poly- 
mer. The effect of the Coulomb interactions is then to 
renormalize the effective thickness of the polymer |17[ I8| . 
However synthetic polymers and single stranded DNA 
both have an intrinsic persistence length of the order 
ip ~ Inm |l9| which could be small compared to the elec- 
trostatic screening lengths. In this paper we explore the 
infiuence of topological constraints on charged polymers 
in cases where the screening length is large or compara- 
ble to the intrinsic persistence length. In Sec. II we start 
by considering the idealized case of unscreened Coulomb 
interactions. This case demonstrates that under long- 
range interactions the topological constraints are pulled 
into tight knots. As discussed in Sec. Ill, this conclu- 



* Electronic address: 
"I" Electronic address: 



panl.dommersnes@phys.ntnu. no 



iantor@post.tau.ac.i 



sion has to be re-examined in real systems due to finite 
rigidity of the polymer, thermal fluctuations, and, most 
importantly, finite screening. Surprisingly, we find that 
tight knots are rather resilient: They remain as global 
equilibrium solutions as long as the overall shape of the 
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polymer is dominated by the (screened) Coulomb inter- 
actions. Tight knots can also remain as metastable states 
for shorter screening lengths, as long as the electrostatic 
bending rigidity is larger than the intrinsic one. Such 
long-lived tight knots have strong influence on the relax- 
ation dynamics of the polymers as discussed in Sec. IV. 
In particular, we find that the most likely way for elimi- 
nating topological entanglements is by diffusion of tight 
knots along the chain; interestingly stronger Coulomb in- 
teractions lead to tighter knots that are less mobile. 



II. UNSCREENED INTERACTIONS 





FIG. 1: The initial (left) and equilibrium (right) conforma- 
tions of a 64- monomer charged polymer, at T = 1.4, forming 
a trefoil knot. (The right figure is reduced by a factor of 2). 



We first consider a simple model of a charged poly- 
mer in which monomers repel each other via unscreened 
Coulomb interactions. The interaction between two 
charges e, in a solvent with dielectric constant £, sep- 
arated by distance r is e^/er, and consequently the over- 
all electrostatic energy of a polymer of N monomers is 
Vc = (e^/e) Xli>j ~ ^jl^ where is the position of 
i-ih monomer. Given a typical separation between ad- 
jacent monomers of a, it is convenient to introduce the 
energy scale Cq = e'^/ea. Initially, we focus on configu- 
rations in which the monomers are locally stretched to 
form smooth straight segments, gradually curving at a 
larger length scale R set by the overall shape. For such 
configurations the Coulomb energy has the form 
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where c is a numerical constant of order unity, and we 
note the following: 

• For any smooth curve, the integral of the 1/r- 
potential leads to a logarithmic divergence, and 
consequently the energy of the polymer is overex- 
tensive, and, consequently, the tension on the poly- 
mer increases as \nN. Therefore, thermal fiuct na- 
tions are irrelevant for a sufficiently long polymer, 
whose shape is determined by minimizing the en- 
ergy. 

• The second term in Eq. ^ can be regarded as the 
Coulomb interaction between charges (or order N) 
on remote parts of the polymer (distances of order 
R). Since typically R ex Na^ the partition of the 
energy between the two parts is not precise, and 
can be changed by redefining R. 

• The Coulomb interaction prefers to keep the 
charges far apart, and the polymer minimizes its 
energy by assuming a shape with maximal R. Thus, 
open polymers simply form straight lines, while un- 
knotted ring-polymers form circles. 

The above argument can be misleading in the case of a 
knotted polymer, as illustrated in Fig. [l[ Here, we used 
Monte Carlo (MC) off-lattice simulations to determine 



the shape of knotted polymers at finite temperature. Our 
model polymer consisted of hard sphere monomers con- 
nected by "tethers" that have no energy but limit the 
distance of a connected pair to 1.05 of the hard sphere 
diameter a. Fig. depicts the results of a simulation for a 
trefoil knot: As an initial conformation (left) in this sim- 
ulation (as well as in the subsequent simulations of more 
complex knots) we used a harmonic representation |2^ 
in which coordinates of the monomers are given as poly- 
nomials in cos(t) and sin(t), where t parametrizes the 
curve. (This provides a relatively clear visualization of 
the knot.) Since the hard core and tether potentials do 
not have an energy scale, the temperature T appears in 
the simulations in the combination ksT /cq^ which we will 
denote as dimensionless temperature T. All simulations 
described in this section were performed for T = 1.4. 
It is customary to represent the strength of the electro- 
static potential by the Bjerrum length is = e'^/ekBT. 
(In water at room temperature -^B=0.7nm.) In our nota- 
tion, the Bjerrum length is simply related to the dimen- 
sionless temperature by Ib = a/T. (Note that for the 
moderate values of A/" = 64 used in this simulation, the 
polymer shape on the right of Fig. [l] is somewhat 'wig- 
gly'; an effect that should disappear for ^ oo due to 
the overextensivity of the energy.) 

Figure |l| clearly shows that in equilibrium the trefoil 
assumes an almost circular shape, with the topological 
details concentrated on a very small portion. (The scale 
of the right side part of Fig. || has in fact been reduced 
by a factor 2 relative to the left figure, and the actual 
linear extent of the equilibrated knot is almost twice its 
initial size.) This behavior can be explained by com- 
paring the long and short-ranged contributions to the 
Coulomb interaction: By expanding its radius, the long- 
range part of the Coulomb energy is reduced by a factor 
of S{N'^/R) oc N. This comes at the cost of bringing 
several charges close together in the tight portion, but 
the latter energy is independent of TV, and can be easily 
tolerated for sufficiently long polymers. 

Because of the highly curved portion, Eq. |l| does not 
apply to tight knots. For a semi-quantitative under- 
standing of the tension that creates such objects, con- 
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FIG. 2: A closed loop (A/'-monomer polymer) folded into 
a shape that can be approximately described as two circles 
consisting oi N — n and n monomers, separated from each 
other by a distance of order n monomer sizes. 



sider a simpler example of an A/'-monomer closed chain 
folded into a shape consisting of a large circle of N — n 
monomers, and a small loop of n monomers, as depicted 
in Fig. ||. For n <C A^, the electrostatic energy can be de- 
composed as Ec{N — n) -\- Ec{n) -f- Ei, where Ec is given 
in Eq. |l|, while Ei is the interaction energy between the 
small loop and the large circle. Assuming that the curved 
strands are separated by a distance of the order na, the 
latter is of the order of 2eon[\n{N/n) + c'], where c' is 
a constant depending on the details of the shape. The 
leading n-dependent part of the total energy is then 



E{N,n)-eon\n\^ — j , (2) 

representing a tension that grows logarithmically with 
the length of the polymer. This conclusion is not lim- 
ited to the shape depicted in Fig. ^ but should apply 
to any smooth linear curve consisting of two portions of 
very different sizes. Equation || thus indicates that from 
purely electrostatic energy considerations n should take 
the smallest possible value, as indeed happens in the case 
of a tight knot in Fig. 0. 

The tightness observed for a trefoil knot also occurs 
in more complicated topologies. Figure § depicts the re- 
sults of equilibration of 128-monomer polymers beginning 
from a harmonic shape on the left, to equilibrium shapes 
(on the right) at T = 1.4. Below each figure we indicate 
the type of the knot in the standard notation C/c, where C 
is the minimal number of crossings the knot can have in a 
planar projection |2^. Since for a given number of cross- 
ings there can exist several different knots, an additional 
subscript k labels the standard ordering of these knots. 
(For C = 3 and 4 there is only one knot, while for C = 8 
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FIG. 3: The initial (left) and equilibrium (right) shapes of 
knots formed by 128-monomer polymers at T = 1.4 (Ib = 
0.7a). A selection of prime knots of varying degrees of com- 
plexity is depicted. (The figures in the right column have 
been scaled down.) The numbers in the left column are the 
standard notations for knot types. 
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FIG. 4: "Coulomb factorization" of composite knots on a 128- 
monomer polymer at T = 1.4. Original (left) and equilibrium 
(right) configurations (scaled down) are shown. 



there are 21 distinct knots [^.) Despite the increasing 
topological complexity of the knots in Fig. |[ their even- 
tual (collapsed-knot) state is reliably represented by the 
semi-quantitative description based on the energetics of 
Fig. |. 

The above arguments indicate the energetic advantage 
of compressing any indivisible topological constraint into 
a tight shape (as opposed to leaving it as an expanded 
structure). However, similar considerations suggest that, 
if possible^ any concentrated region of charge should split 
into smaller elements placed as far as possible from each 
other. Such a reduction is not possible for the prime 
knots considered in Fig. |3[ which (by definition) cannot 
be separated into several parts connected by a single line. 
In contrast, composite knots are formed by joining several 
prime factors together, and Fig. || presents initial and 
final (equilibrium) states of several such knots on 128- 
monomer polymers. The notation below each knot indi- 
cates its constituent prime components. The Coulomb 
interaction clearly "factorizes" any composite knot, sep- 
arating its elements as far as possible. However, since the 
typical interaction energies between the prime factors are 
only a few eo, thermal fluctuations (T = 1.4) in the dis- 



tances between these tight regions are quite pronounced. 



III. BEYOND 'IDEAL' KNOTS 

Many of the results in the previous section are in fact 
known to knot theorists, who have investigated long- 
range repulsive interactions with the aim of finding a 
knot-invariant energy |2^, The basic question is 

whether a properly scaled energy of the ground state con- 
figuration (the ideal state) for certain choices of interac- 
tion functions can be used as a means of distinguishing 
different knot types. An example of such an interaction 
is Simon's 'minimal distance' between the strands func- 
tion, or a repulsive 1/r^ type interaction which pro- 
duces symmetric spread out ground states. In Ref. it 
was conjectured that minimizing knot- invariant energies 
should decompose a knot into prime sub-knots and simu- 
lations with 1/r^ interactions support this Electro- 
static interactions do not generate useful knot- invariant 
energies, since, in the absence of excluded volume in- 
teractions, knots on a continuous curve are collapsed 
to a point providing no (cut-off independent) way 
of identifying knots. (Indeed, in the simulations of the 
previous section knots were tightened into compact ob- 
jects whose extent was determined by the monomer size.) 
While this conclusion may be disappointing to a knot- 
theorist, it is encouraging from the perspective of poly- 
mer science, since it is easier to describe the properties 
of tight entanglements, without having to worry about 
their precise topology. However, this is the case only if 
we can demonstrate that tight knots survive for realistic 
polymers subject to electrostatic interactions in actual 
solvents. Accordingly, in this section we shall include 
additional attributes present in such situations, and con- 
sider the effects of bending rigidity, thermal fluctuations, 
and (most importantly) of a finite screening length. In 
these circumstances the size of the knot can be signifi- 
cantly larger than in its maximally tight state; neverthe- 
less, tight knots can still remain. 



A. Bending rigidity 

Many microscopic aspects of polymers are captured at 
a mesoscopic scale by a curvature energy, describing its 
resistance to bending. In a charged polymer one should 
distinguish between the intrinsic bending rigidity, and 
an effective rigidity which includes the electrostatic con- 
tributions. The latter arises because bending a straight 
segment brings the monomers closer and thus increases 
the Coulomb energy. The former can be represented by 
a length ip at which, in the absence of other interac- 
tions, the transverse thermal fluctuations of the poly- 
mer become of the same order as the length-scale itself, 
or at which orientations of the bonds become uncorre- 
lated. Simple analysis relates £p to the bending rigidity 

and temperature by = ksTip. In charged polymers. 
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Rij ^ aN^ with u ^ 0.59. Consequently, the linear ex- 
tent of the entire polymer is approximately Ri){N/Ni)). 
If a weak force / is applied to a segment of spatial 
extent i?5, that segment is stretched y by an amount 
X ^ Rlf /ksT. The size of a blob is determined by re- 
quirement that X ^ leading to R^ ~ ksT/ f. An 
open charged polymer can also be viewed as a stretched 
chain formed from such blobs [|9[ |3^, while a ring- 
polymer is a circle of such blobs. The force stretching 
a blob in an object of this type is €oa{Ni)/ Ri^)^ ln(A^/A^b). 
By substituting this force into the expression for blob 
size, and solving it, we extract the number of monomers 
in each blob as 



FIG. 5: Equilibrium configurations of a 128-monomer trefoil 
knot: (a) a tight 20 monomer) knot at T = 5, and (b) an 
expanded 60 monomer) knot at T — 10. 



£p should be measured in the presence of high salt con- 
tent, so that electrostatic contributions to rigidity are 
screened out. It is reasonable that the bending rigid- 
ity, rather than monomer size, should determine the size 
of a tight knot. The energy for bending a segment of 
length ^, with radius of curvature also of order of is f<i/i 
with a dimensionless shape-dependent prefactor. For the 
shape depicted in Fig. ||, there is now a bending cost of 
Eb ^ nina which competes with the electrostatic energy 
in Eq. g. By minimizing the sum of these energies, we 
find that the optimal knot size is 



eoHN/rik) y CoHN^Co/^) 



(3) 



where we have omitted numerical prefactors of order 
unity. This result indicates that the knot in stiff polymers 
of moderate size N can be as large as hc/cq = \/^p/^B^ 
and becomes compact only for N ~ ex.-p{f<i/eo). 



B. Thermal Fluctuations 

At high temperatures, entropic factors (which fa- 
vor crumpled states) compete with electrostatic effects. 
While the latter dominate on sufficiently long length- 
scales, at short length-scales fluctuations are important. 
This competition can be visualized by a simple blob pic- 
ture [|l|. If a strong external force / is applied to a self- 
avoiding polymer without electrostatic interactions, it is 
stretched to a linear form. This linear object, however, 
has a finite width i?5, and can be regarded as a chain 
of blobs of this size. On length-scales shorter than the 
blob size, the external force has negligible effect, and we 
can relate R^ to the number of monomers forming the 
blob via the usual relation for self-avoiding polymers [ffl: 



eoHN/N^) 
ln7V(a/£B)^ 



ln(7V/f — ) 



(4) 



Of course, the blob picture is meaningful only if A/'^ is 
larger than unity. Thus blobs can appear only for tem- 
peratures T InA'; and for A^ = 128 we expect to see 
the blobs for T ^ 5. Fig. || depicts equilibrium shapes 
of a trefoil knot at T = 5 and T = 10, and we see the 
appearance of a wiggly structure in the higher tempera- 
ture regime. At such high temperatures, we expect knots 
to have a size typical of that in a non-charged polymer 
consisting of A^& monomers. The exact size of the knot re- 
gion in non-charged polymers in three-dimensional space 
is not known; simulations suggest that knots are local- 



ized ll2 



^ L3|, but not compact [y^. The size of the blob 
in Fig. |5| is too small for any kind of quantitative study, 
but we clearly see that the knot is no longer maximally 
compact. 



C. Screened Interactions 

A charged polymer in solution is accompanied by neu- 
tralizing counter ions, and potentially other charged ions 
due to added salt. In general, the effect of these ad- 
ditional ions on the charged polymer is quite compli- 
cated, and dependent on the intrinsic stiffness, strength 
of the charge, and valency of counter ions However, 
in many cases the net effect can be approximated by 
a screened Coulomb potential V = (e^/gr) exp (— r/A), 
where A is the Debye screening length l32| . Since the pre- 
vious arguments for the tightness of charged knots rely 
on the long-ranged part of the Coulomb interaction, we 
may well question if and when tight knots survive with 
screened forces. 

It is important to realize that Coulomb interactions 
affect the polymer on scales much larger than A, due 
to increased bending rigidity. Curving a straight poly- 
mer to a radius R brings its charges closer, resulting 
in an extra energy cost of order {e^ /eR)}? for screened 
Coulomb interactions, where A = A/a is the reduced 
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FIG. 6: Radius of gyration, i?g, of the ground state configu- 
ration of a trefoil knot, as a function of the screening length 
A for a 128-monomer polymer. Rq has been normalized by 
the length of the polymer Na. 



screening length. This can be regarded as an effective 
bending rigidity, which (in the presence of thermal fluctu- 
ations) leads to the Odijk-Skolnick-Fixman persistence 
length H of 4 = X^e'^/ieksTa'^) = ^bA^. The electro- 
static persistence length is in general much larger than 
the screening length. In terms of our reduced variables 
£c = ^c/a = X'^ /T. This expression is valid provided that 
the length scales considered are larger than the screening 
length, and T < A. 

For very large A, comparable to the size of the poly- 
mer, the effects of screening are not very important: E.g., 
Eq. H for electrostatic energy of a knot remains valid if 

is replaced by A, and similar replacements should be 
made in Eq. |^ for the knot size in a stiff polymer, or 
m Eq. I for the blob size. In all these expressions, the 
number of monomers enters only in a logarithm, and, 
consequently, its replacement by A does not significantly 
change the result. Eq.J| for the optimal knot size is valid 
(with N replaced by A) only if the knot is smaller than 
the screening length. This condition, A > an/c, leads to 
the cross-over boundary 




(5) 



which is equivalent to ic > ip. We thus conclude that 
a tight knot can exist only when the overall bending 
rigidity is dominated by electrostatic contributions. For 
smaller values of A, the short range repulsion can no 
longer bend the knot into a tight shape. 

Note that the analysis leading to Eq. ^ only demon- 
strates the local stability of a tight knot. The global 
energy minimum could still occur for a spread out con- 
figuration. To decide on the latter requires estimates of 
the energy difference between the two configurations, and 
depends on microscopic details, as well as the length of 



the polymer. A circle with a tight knot, and the spread 
out knotted shape, both have a bending energy (at large 
scales) of the order of kBTic/R- Since the circular shape 
has a larger radius, it has a lower energy, the energy 
difference scaling as k^Tl^/ {Na)^ if both radii are pro- 
portional to the polymer length. The tight knot in the 
former has an additional local energy cost, which is of 
the order of kBT^l^/ a) (possibly with logarithmic cor- 
rections), but independent of N. Thus, we expect the 
configuration with a spread-out knot to have a lower en- 
ergy only for < ^B, i-e. for screening lengths 



A < Ac ^ aViV. 



(6) 



Note that the limiting value of A still corresponds to a 
persistence length of the order of the extended polymer, 
i.e. the polymer shape is determined by energy consid- 
erations, and thermal fluctuations have little effect, at 
this point. We verified this conclusion by numerically 
determining the shape of the trefoil that minimizes the 
screened Coulomb interactions. Figure || shows the ra- 
dius of gyration as a function of the screening length. For 
screening lengths larger than Ac ~ ^AaN^^^^ the knot 
switches from a loose to a tight configuration. 

Let us briefly explore the possibility of tight knots in 
nucleic acids. Double stranded DNA has a bare persis- 
tence length of rsj 50nm, which is much larger than 
typical screening lengths, and consequently is not likely 
to incorporate any knots tightened by Coulomb interac- 
tions. However, measurements on single stranded DNA 
in high salt concentrations [|l9| suggest a much smaller 
intrinsic £p ~ Inm, and presumably a similar (or smaller) 
value applies to single stranded RNA. Tight knots should 
then occur for single stranded nucleic acids for reasonable 
screening lengths of the order A ~ lOnm. This could 
for example be relevant to the experiments of Ref. [Q, 
where artificial knots in single stranded RNA were used 
to demonstrate the existence of a topology changing en- 
zyme. Knotted polymers are often distinguished from un- 
knotted ones by electrophoresis However, if the knot 
is tight, the knotted polymer may have an electrophoretic 
mobility close to that of a ring polymer, making such de- 
tection problematic. 



IV. TIGHT KNOTS AND DYNAMICS 

Tight knots are created whenever a polymer is under 
tension; the source of tension need not be long-range re- 
pulsions. For example, it has been argued that tight 
molecular knots appear in polymer systems undergoing 
crystallization, as crystallization at one point may create 
tension in other parts of the chain [Q. Polymers in a 
strong shear flow are also subject to tension [^, and 
may even undergo a coil-stretch transition as a result [|6| . 
It is plausible that stretching could tighten loose knots 
in the chain. Once created, such molecular knots should 
be quite stable and thus account for long-time memory 



effects observed in polymer melts Q. However, molecu- 
lar dynamics simulations suggest that once the tension is 
removed a tight knot opens up in a short time . With- 
out being systematic, here we examine a couple of dy- 
namical issues pertaining to charged tight knots, namely 
their creation in a high temperature quench, and their 
relaxation by diffusion along the chain. 



A. Tightening by Quench 

It is quite likely that when topological entanglements 
are first formed, e.g. in the process of cyclization of a 
polymer, they are spread out over the whole chain. Sub- 
sequent tightening then occurs upon increasing tension. 
In the case of charged knots, this process is illustrated in 
Fig. 0. Here, the initial configurations are the spread out 
harmonic representations, which soon evolve into loops 
separated by tight elements. The relaxation process then 
slows down as one of the loops grows at the expense of 
the others. A universal last stage is the appearance of a 
structure reminiscent of Fig. ||, with two loops separated 
by a tight 'slip-link.' We observed the same sequence in 
simulations where the initial configurations was an equi- 
librated (random walk) knot. The formation of the two 
loops separated by a slip-link was again relatively fast, 
and the rate limiting step was the sliding of one loop 
through the tightly packed monomers at the slip-link. 

B. Diffusion of Tight Knots 

As demonstrated in the previous situation, tight knots 
slow down the relaxation of the polymer to its eventual 
equilibrium state. Here we study such relaxation more 
explicitly for a knot in an open charged polymer. In this 
case there is no topological constraint, and the polymer 
is expected to unknot to achieve its equilibrium state. 
Does a tight knot in an open chain relax by becoming 
loose and opening up, or by sliding (diffusing) to one end. 
As demonstrated in Fig. ||, the latter is the case: The 
initial configuration (in a chain of = 64 monomers with 
unscreened interaction) remains tight, indicating that the 
stretching force from the monomers at the ends of the 
chain is larger than from those forming the knot. In the 
simulation, the knot's position fiuctuates for some time in 
the middle, before moving to one direction. The eventual 
unknotting occurs when the diffusing tight knot reaches 
the end of the polymer. 

A tight knot in the middle of an open chain is in 
a meta-stable state. We can estimate a potential en- 
ergy for the tight knot by considering a charge Q = ne 
along a charged chain of N monomers. The Coulomb 
energy then depends on the position of this charge A^i, 
as E = kBT{£B/a)nln[Ni{N - Ni)]. This energy is min- 
imal when the charge Q is at either endpoint of the poly- 
mer, i.e. for A/"! = or N. Note that the force pushing 
the extra charge towards the end scales with -^b, and we 
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a) b) 



FIG. 7: Time evolution (using Monte Carlo dynamics) of (a) 
3i and (b) 819 knots, from the initial (harmonic) geometry 
(top) through an intermediate state when the knot "strangles" 
the loop close its middle, and to a final state (bottom) when 
the knot is localized. A similar sequence takes place for all 
other prime knots in the simulations of Fig. H. 



may naively expect that the resulting relaxation becomes 
faster as the Coulomb energy is increased. In fact, the 
opposite occurs for charge knots, with relaxation slowing 
down as Coulomb interactions become more dominant. 
The reason is that increased charging energy leads to 
a higher tension and more closely packed monomers in 
the knot. Any motion of the knot requires some inter- 
nal rearrangements of these monomers, accompanied by 
pulling in some monomers from the straight portions of 
the chain. This necessitates overcoming an energy bar- 
rier of ~ -^B In A^, and consequently higher charged knots 




FIG. 8: Unknotting of a charged polymer with N — 64 
monomers and unscreened Coulomb interaction of strength 
T — lA (Ib — 0.7a). The initial configuration is a tight knot 
in the middle of the chain. Rather than open up gradually, 
the knot slides along the polymer and remains localized until 
it reaches the end. 



are tight and harder to move. Since rearrangements re- 
quire a large activation energy, the knot remains stuck in 
position. This is quite similar to what happens to a knot 
in a polymer under strong tension [|4|. 

While with unscreened Coulomb interactions the tight 
knot feels a potential that drives it to one end, there is no 
such force when the interactions are screened (unless the 
distance between the knot and the endpoint of the poly- 
mer is of the order of the screening length). The energy 
barrier preventing the loosening of the knot is also finite 
in this case. The resulting dynamics for a chain of 128 
monomers with a screening length of A = 6a is demon- 
strated in Fig. ^; despite the screening the knot remains 
tight until it diffuses to one end. The characteristic time 
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X = 6a 




FIG. 9: Monte Carlo dynamics of a tight knot in a chain with 
iV = 128 monomers and screened interactions. The screening 
length is A = 6a, roughly the size of the knot. The knot shows 
no sign of opening up, it remains tight till it reaches the end 
of the polymer. 



scales for the relaxation of the knot can be estimated 
as follows. The time for diffusion over a distance Na 
scales as a^N'^ / Dknot^ with the knot diffusion coefficient 
behaving as Dknot oc I) exp (— f^/ZcBT). Here, D is the 
diffusion constant for a single monomer, while the activa- 
tion energy for local rearrangements necessary for motion 
of the tight region is roughly Ed ^ k^T{i^/a) log(A/a). 
There is also the possibility that the knot becomes loose, 
escaping the local minimum of the tight configuration. 
The energy barrier for the latter is « /c^T-^eA/a^, with 
a corresponding time scale of r ~ {a^ / D) ex.-p {Sb / kBT) . 
In time r, the knot can diffuse a distance L ^ ^/Dknot^- 
We thus estimate a "processivity length" over which a 
tight knot diffuses, before opening up, by 

Lp (X aexp (C ^sA/a^), (7) 

where C is a constant of order unity. The processivity 
length increases strongly with the screening length A and 
quickly reaches a macroscopic length, indicating that the 
relaxation of a tight knot will be by diffusion along the 
chain, even for very long chains. Also note that Lp is 
in general much larger than the electrostatic persistence 
length which only grows quadratically with the screening 
length (4 ~^BA/a2). 

V. DISCUSSION 

We have shown that long-ranged Coulomb forces gen- 
erate a tension that tightens topological constraints into 
dense localized regions, leaving the rest of the polymer 
unentangled. For knots on ring polymers, we confirm the 
"factorization" of composite knots into their prime com- 
ponents. Tight knots remain, even when the Coulomb in- 
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teraction is screened, as long as the electrostatic contribu- 
tions dominate the rigidity of the polymer. Once formed, 
tight knots drastically slow down the equilibration of the 
polymer (or polymer solution), as they typically relax by 
diffusion along the backbone. If the Coulomb interac- 
tions are strong enough, the knot is pulled so tight that 
it is unable to diffuse, and its position appears frozen. 
This is different from uncharged polymers where molecu- 
lar dynamics simulations in ref. ||3^ find that tight knots 
in short uncharged polymers open up rapidly. Our re- 
sults predict that tight knots in polyelectrolytes can be 
very stable and cause long relaxation times. While we 
have focused on single polymers, it is natural to specu- 
late about similar behavior in solutions of many chains. 
It is indeed quite likely that inter- chain entanglements 
are also tightened in polyelectrolyte solutions and gels. 

Additional consequences of tight knots are in their 
influence on mobility (electrophoresis), and on the me- 
chanical strength of polymers. It has been shown re- 
cently by direct measurement on DNA and actin fil- 



aments that knots significantly weaken the strand 
Similarly, molecular dynamics simulations of knotted 
polyethylene chains also find that the strands becomes 
weaker, and typically break at the entrance point where 
the straight segment ends and the tight knot begins . 
Single stranded DNA is relatively fragile and sometimes 
breaks during electrophoresis or when subject to flow; 
tight knots may well be responsible for this phenomenon. 
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